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Abstract 
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We introduce a new variational characterization of Gaussian diffusion processes 
as minimum uncertainty states. We then define a variational method constrained by 
kinematics of diffusions and Schrodinger dynamics to seek states of local minimum 
Q_i' uncertainty for general non-harmonic potentials. 

> 
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■ 1. Introduction 

Nelson stochastic quantization, originally proposed in 1966 0, is currently recognized 
as an independent and self-consistent formulation of quantum mechanics in the language 
of stochastic processes. 

The original Nelson scheme, containing some ad hoc assumptions, has been later re- 
considered and conceptually clarified in the classic work by Guerra and Morato ||. They 
proved that quantum dynamics can be obtained by a stochastic variational principle, i.e. 
by suitably extremizing the classical action along diffusive trajectories replacing the deter- 
ministic ones. In this way stochastic variational methods have creeped in into quantum 
physics. Both the practical and the conceptual advantages of the variational scheme have 
since been exploited in a number of different contexts 0, ||, 0. 

In this work we introduce a variational method that yields a new characterization of 
Gaussian diffusion processes. The analysis is carried out for the one- dimensional case, but 
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the results can be easily extended in any number of dimensions. We prove that a necessary 
and sufficient condition for a diffusion process q(t) to be Gaussian is that a certain positive- 
defined functional of the probability density p(x, t) acquires a non-zero global minimum 
at each fixed instant of time. This functional is built by considering the variances Ag 2 of 
q(t) and Au 2 of the logarithmic derivative u(x,t) = ud x \\.np(x,t)], where v is the diffusion 
coefficient. We then consider the "uncertainty" product Aq 2 Av 2 as the quantity to be 
extremized against smooth variations 5p of the probability density. 

For the Nelson diffusions of stochastic mechanics the quantity u(x,t) is the so-called 
osmotic velocity, and Aq 2 Au 2 is the osmotic uncertainty product. One can then rephrase 
our result by saying that Nelson Gaussian diffusions are all and only those that minimize 
the osmotic uncertainty product at each fixed instant of time. 

Historically, an uncertainty principle involving the diffusion coefficient had been already 
conjectured by Furth in the study of Brownian Motion ||, and explicitly derived for 
the Nelson diffusions of stochastic mechanics by exploiting Schwartz's inequality 0, ||. 
Saturation of the osmotic uncertainty then yields the standard harmonic-oscillator coherent 
and squeezed states || . Obviously, the same states can also be recovered in the stochastic 
variational approach. 

In fact, they are identified in two steps. Knowledge of the Gaussian density allows 
first to derive their phase. This is achieved by imposing Fokker-Planck kinematics through 
the continuity equation. The latter connects the density with the gradient of the phase, 
i.e. the so-called current velocity v. The harmonic-oscillator potential is then determined 
by means of Schrodinger dynamics in the Hamilton- Jacobi-Madelung hydrodynamic form. 
Thus, in this scheme, the kinematical and dynamical inputs are exploited a posteriori, once 
the osmotic uncertainty functional has been minimized. 

On the other hand, in stochastic mechanics there is a current uncertainty term Aq 2 Av 2 
that has to be added to the osmotic one to yield the full Heisenberg uncertainty product 
Aq 2 Ap 2 . Minimization of the osmotic term alone takes care of the current uncertainty by 
fixing it to be either zero (Heisenberg minimum) or related to the time-variation of the 
wave packet spreading Aq (Schrodinger minimum). 

It is however possible to consider a variational extremization for the sum of both parts, 
which is in fact a functional of the density and the current velocity. The latter are not 
independent but are related by Fokker-Planck kinematics and Schrodinger dynamics, so 
that minimization can be correctly performed only by including them a priori as constraints 
in the variational procedure. Furthermore, the osmotic uncertainty product Aq 2 Au 2 is a 
quantity defined at a fixed instant of time. Upon integration over an arbitrary time- 
interval [ti,^]) if is mapped into a dynamical functional of the process (intuitively, a 
sort of uncertainty "action"). The correct procedure to minimize this dynamical osmotic 
uncertainty is then again to impose the kinematical and dynamical constraints. 

In the present work we first develop the free variational minimization of the osmotic 
uncertainty, a result that holds for any diffusion process. We then carry out the study for 
the various constrained variational schemes. 

It is clear that the constrained approach to the minimum osmotic uncertainty, by re- 
stricting the class of allowed variations, should in principle yield a larger class of solutions. 
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In particular, it should allow to assess the existence and to determine the relative minima of 
the Heisenberg uncertainty product as functions of the assigned external potentials. Unfor- 
tunately, due to the intrinsic mathematical complications, we are at present able to carry 
out the analysis and present a complete solution only for the kinematically constrained 
problem. 

2. Variational characterization of Gaussian diffusion processes 

Consider a one-dimensional diffusion process q(t) with the associated normalized proba- 
bility density p(x, t) and a diffusion coefficient v which may be constant or time-dependent. 
Introduce the quantity 

u(x,t) = ud x [\np(x,t)\ = (1) 

and consider the following functional of p(x,t): 

F[p] = Aq 2 Au 2 , (2) 

where Aq 2 = Var{q} = (q 2 ) — (q) 2 , and Au 2 = Var{u] = (u 2 ) — (u) 2 are the variances 
of q(t) and u(x,t) as functions of time. Here and in the following ( • ) denotes expectation 
with respect to the density p(x,t): for u(x,t) as function of the process, 

/+oo 
d x p(x,t)dx = 0, (3) 

since the density vanishes at infinity. Thus u(x,t) always has, by construction, zero expec- 
tation value, and Au 2 is simply (u 2 ). 

We observe that F is a smooth, positive-defined functional of the density p. Pictorially, 
Aq 2 and Au 2 are in competition: the first quantity measures the width of the probability 
density, while the second measures its sharpeness; we therefore expect on intuitive grounds 
that there exists an absolute minimum of F. To verify this assertion and to determine the 
explicit minimum value of the functional we extremize F under smooth variations 5p of 
the probability density; they are taken to be vanishing at the boundaries of integration. 

Since minimization of F must preserve normalization of probability, the variations will 
be constrained by the requirement that at each fixed instant of time / p(x,t)dx = 1. We 
thus consider as the functional to be varied 

/ + OC 
p(x,t)dx. (4) 

The Lagrange multiplier \(t) shall be determined a posterioribj taking the expectation of 
the variational equation. The latter is obtained by imposing the vanishing of the variation 
SJ 7 of the uncertainty functional: 

/+oo 
5p(x,t)dx = 0. (5) 
-oo 
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Long but straightforward calculations lead to the following variational equation for 

u(x, t): 

Au 2 (x 2 - 2(q)x) - [u 2 {x, t) + 2ud x u(x, t)] Aq 2 + A(t) = . (6) 
Taking the expectation of both sides one finds 

X(t) = Au 2 [(q) 2 - 2Aq 2 ] . (7) 

Replacing the above expression for A in eq.(6) gives the following inhomogeneous nonlinear 
differential equation of the Riccati type for u(x, t): 

Aq 2 [u 2 (x, t) + 2ud x u(x, t)] = Au 2 x 2 - 2Au 2 (q)x + Au 2 [(g) 2 - 2Aq 2 ] . (8) 

The right hand-side is of the form a(t)x 2 + f3(t)x + "f(t): it follows that the solution 
u(x,t), if it exists, must be of the form a(t)x + b(t). The condition that u(x,t) must have 
a vanishing expectation fixes b(t) = —a(t){q), so that the solution is finally of the form 

u(x,t) = a(t)[x-(q)}. (9) 

The coefficient a(t) is determined by the requirement that u satisfies identically the 
variational equation; inserting expression (9) in eq.(8) and equating, respectively, the co- 
efficients of the powers x 2 ,x and x°, yields a(t) = —Au/Aq and 

Aq 2 Au 2 = v 2 . (10) 

From the explicit form 

«M) = -(?)) > (11) 

we readily trace back the expression for the normalized probability density 

**■" = w^M JJ ^f} ■ (12) 

It is possible to check, by convexity of the functional or by taking second variations, 
that the extremum (10) is indeed a minimum. 

In conclusion, we have proven the following 
Theorem: A necessary and sufficient condition for a smooth probability density p(x, t) 
to be Gaussian is that it makes stationary the functional Var{q} ■ Var{u}, where u = 
ud x [\n p(x,t)} and v is an arbitrary function of time (for diffusion processes, the diffusion 
coefficient). The extremum corresponds to a minimum with value v 2 . 

Sufficiency was just shown; necessity trivially follows from the Gaussian characteriza- 
tion, eqs.(ll)-(12), and from the definition of F[p\. 
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3. Minimum uncertainty Nelson diffusions: coherent states 

The theorem proved in the previous section finds application in the framework of 
stochastic formulation of quantum mechanics. We recall that Nelson stochastic quanti- 
zation []TD associates to each quantum state of a point particle of mass m a configurational 
diffusion process q(t) governed by Ito's stochastic differential equation 



dq(t) = vt + )(q(t),t)dt+ dw(t) , dt > . 

\ 2m ) 



(13) 



Here v^(q(t),t) is the forward drift, h is Planck constant, v = h/2m is the diffusion 
coefficient, and dw(t) is the time-increment of a Wiener process w(t), with expectation 
(dw(t)) = and covariance (dw 2 (t)) = 2dt. The forward and backward drifts, vr + )(x, t) 
and f(_)(x, t), are defined by the conditional expectations 



v {+) (x,t) 



lim 



q(t + At) - q(t) 



vi\(x,t) = lim 



At 



q(t) - q(t - At) 



q(t) = x) , 



(14) 



At 



q(t) = x) . 



They represent respectively the mean forward (backward) velocity fields. 

In the hydrodynamic picture of the process, the drifts are replaced by the osmotic 
velocity u(x, t), 



u(x, t) 



— d x [ln p{x,t)} . 



and by the current velocity v(x,t), 

v(x, t) 



V M (x,t) +V(-)(x,t) 



(15) 



(16) 



Finally, Fokker-Planck equation for the probability density p(x, t) takes the form of the 
continuity equation 

d t p(x,t) = -d x [p(x,t)v(x,t)) . (17) 

At the dynamical level to each single-particle quantum state ^f(x, t) written in the 
hydrodynamic form 

ty(x,t) = 



p(x, t) exp 

there corresponds in stochastic mechanics the diffusion process q(t) with 



-S(x,t) 



and 



p(x,t) = \*{x,t)\' 



v(x,t) = — d x S(x,t) . 
m 



(19) 
(20) 



The complex Schrodinger equation with potential V(x,t) for the wave function \l/ is 
then equivalent to two coupled real equations for the probability density p and for the phase 
S (or, alternatively, for the osmotic and current velocities u and v) . They are, respectively, 
the continuity equation (17) and the gradient of the Hamilton- Jacobi-Madelung (HJM) 
equation which, in terms of u and v , reads 

md t v + mvd x v — mud x u — —d 2 u = —d x V(x, t) . (21) 

The correspondence between expectations and uncertainties defined in the stochastic 
and in the canonic formulations of quantum mechanics are || 

<?)* = (<?), (p)« = m(v), 



Aq = Aq, (22) 

Ap 2 = m 2 [Au 2 + Av 2 ] , 

where q and p are the position and momentum operators in the Schrodinger picture, ( ■ )<t 
denote the expectations of the operators in the given state *ff, ( • ) denote the expectations 
of the stochastic variables in the state {p, v} in the Nelson picture, and A(-) denote the 
variances. 

From eq.(22) it follows that in the Nelson picture the Heisenberg uncertainty product 
Aq 2 Ap 2 = H[p,v], a functional of the density p and of the current velocity v: 

H[p,v] = m 2 Aq 2 Au 2 + m 2 Aq 2 Av 2 . (23) 

Apart from the constant factor m 2 , the first term in the right hand side of eq.(23) is 
the osmotic uncertainty functional F[p] introduced in section 2, while the second term 
m 2 Aq 2 Av 2 is the current uncertainty functional. 

We now have the choice to either extremize only the osmotic term F[p] under variations 
Sp, or the total uncertainty functional H[p,v] under variations 5p and 8v. In the former 
case, one obtains the Gaussian structure (10)-(12), and determines the current velocity 
through the continuity equation (17). The result is 

v(x,t) = (v) + [ x -( q )]—JL. (24) 

If we take Aq to be constant in eq.(24), we have v — (v), and Av = 0. The quantum 
uncertainty then collapses to the minimum osmotic uncertainty fr 2 /4, i.e. to the standard 
absolute Heisenberg minimum. The corresponding states have Gaussian density (12) and 
phase exponent of the form (p)x: they are the harmonic-oscillator coherent and squeezed 
states, the potential being determined via the HJM equation (21). 
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With Aq time-dependent we have that Av 2 = (dAq/dt) 2 . It is straightforward to show 
that it is the Schrodinger position-momentum correlation: 



5^=535 [<«*>-»- <»»>]• (25) 

where {qp) symm = (qp + pq)/^- We thus have 

AfAp 2 = + [(gp) symm - (g)(p)] 2 . (26) 

From the above expressions we see that minimum osmotic uncertainty with time- 
dependent Aq yields the Schrodinger minimum of the Heisenberg uncertainty product 
|TT| . For the detailed structure of coherent and squeezed states in stochastic mechanics we 



refer to the work by De Martino et al. 



4. Kinematical constraint and absolute minimum uncertainty 

So far, it has been shown that minimizing the osmotic uncertainty functional F[p] 
without constraints reproduces the whole minimum-uncertainty structure of quantum me- 
chanics provided one exploits, after minimization has been carried out, the kinematics of 
diffusions and the Schrodinger dynamics of quantum states. 

On the other hand, a fundamental aspect of the calculus of variations regards the 
implementation and the meaning of constraints. Minimizing a free functional and then 
requiring compatibility of the obtained results with some constituent equations is in general 
not equivalent to minimizing the same functional with the prescribed equations added as 
constraints. True, it is meaningless to seek a constrained minimization of the osmotic 
uncertainty F[p], since it is a quantity defined and evaluated at each fixed instant of time. 
However, wishing to implement the kinematical and dynamical equations as constraints, 
one has simply to replace F[p] with its time- integral over an arbitrarily fixed interval 
[ti,^]; this new functional is now dependent on the dynamics: it gives the time-evolution 
of the osmotic uncertainty along the process, and it is then sound to seek its constrained 
minimization. 

We begin by showing that the crucial aspects of uncertainty stem from the osmotic 
contribution alone. To this end, we first minimize the total uncertainty functional H[p,v]. 
In this case the need for constrained minimization is imposed ab initio. Namely, p and 
v are linked by the continuity and HJM equations, so that it is impossible to consider 
independent variations 5p and 8v. 

In all generality, one should add both the continuity and the HJM equations to H 
and minimize the resulting constrained functional. However, we start by attacking the 
simpler problem of minimizing H with the kinematical constraint alone. In this way 
dynamics is exploited after minimization, to select the quantum mechanical processes that 
are compatible with the variational results. 

Again, when dealing with time-evolution equations, the quantities to be minimized are 
not the original functionals, but their time-integrals over an arbitrary interval. Keeping as 
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before the normalization constraint, one has to consider the following functional of p(x,t) 
and v(x, t): 

H[p,v] = f t2 dtm 2 Aq 2 (Au 2 + Av 2 )+ f t2 dt / + °° dx\ Xo^p+X^x, t)[d t p+d x (pv)}\ , (27) 

Jt x Jt x J-oo I J 

with the arbitrary functions Xo(t) and Xi(x,t) as Lagrange multipliers. 

By requiring stationarity of TL with respect to independent variations of p and v we are 
left with the following equations: 

2m 2 Aq 2 (v - (v)) - 9 x Ai = , (28) 

m 2 (Au 2 + Av 2 )(x 2 -2(q)x) +m 2 Aq 2 (v 2 -2(v)v - ^^^-) + X (t) -dtX.-vd^ = . 

(29) 

The unknown multipliers Ao and Ai can be eliminated from these relations by inserting 
eq.(28) and its time derivative into the space derivative of eq.(29). Differentiation with 
respect to time obviously yields two distinct results according to whether the spreading 
Aq is taken to be constant or time-dependent. 

In the first case one has, recalling the definition of u, eq.(6), 

d t (v - (v)) + vd x v - Au2 + AV \ x - (q)) + l -d x (u 2 + -d x u) =0. (30) 

Aq 2 \ m J 

For the Nelson diffusion of quantum mechanics the above equation must hold together 
with the HJM dynamics, eq.(21). One can show, after some manipulations, that a sufficient 
condition for the two equations to be consistent is v — (v). This condition reduces eq.(30) 
to the form 

8 X |a«V - 2(g)x) - [u 2 + ^u^jAq 2 + A(t) j = . (31) 

This is just the space derivative of eq.(6) of section 2, with Aq now fixed to be constant, 
and we again recover the Heisenberg minimum uncertainty states. 

In the case of time-dependent Aq, one has instead, from eqs.(28)-(29), 

d t (v-(v)) + 2(v-(v))d t \nAq + vd x v - ^±^! ( x - ( q )) 

Au 

Ti 



+ \d x (u 2 + -d x u\ = 0. (32) 
2 V my 



The problem is now slightly more complicated, but one can convince himself, proceeding 
as before, that HJM dynamics is compatible with the choice of the linear v, eq.(24). By 
inserting this form of v in eq.(32), the linear u of minimum uncertainty is recovered if the 
time-evolution of Aq is ruled by 



One can realize by direct computation that the above expresses the condition for the 
current uncertainty functional Ag 2 Aw 2 , namely the Schrodinger part of the uncertainty 
product, to be a constant of the motion. Eq.(33) can be explicitely solved to give 



Aq(t) = Aq ^l + f3t, (34) 
where Ag = Aq(t = 0), (3 = 2 (dAq/dt) t=0 / Aq(0); consequently, 

h 2 

m 2 Aq 2 Au 2 + m 2 Aq 2 Av 2 = — + m 2 [Aq 2 (dAq/dt) 2 ] t _ Q . (35) 

The corresponding normalized solutions are states of time independent Schrodinger 
minimum. They are characterized by a wave packet width that evolves according to the 
diffusive dispersion law Eq.(34). 

We see that constrained minimization of the total uncertainty functional selects all 
the states of Heisenberg minimum uncertainty, but only the special class of constant 
Schrodinger minima. This is in contrast with the results obtained extremizing the free os- 
motic functional which allowed to derive all the states of minimum quantum uncertainty. 
For this reason, we consider again the osmotic uncertainty and seek for its constrained 
minimization. Proceeding as before, we introduce the functional 



C[p,v] = j 2 dtm 2 Aq 2 Au 2 
Jti 

+ jf dt dx{\ (t)p + \ l (x,t)[d t p + d x (pv)]j . (36) 

Performing the independent variations of p and v and requiring stationarity of C leads 
to a couple of equations that replace eqs.(28)-(29) 

d x Xi = , (37) 

m 2 Au 2 (x 2 - 2(q)x) - Ti 2 Aq 2& ^- + A (t) - d t \ x - vd x \ x = . (38) 

\[P 

Eliminating the unknown multiplier Ai from the equations one is left with 

d x |a«V - 2(q)x) - (u 2 + ^d x v)jAq 2 + A(t) j = . (39) 

We are again in the general case of section 2, eq.(6). We remark that now, contrary 
to the situation faced in the derivation of eqs.(30)-(31), no hypotheses need to be made 
on the time dependence of Aq. In conclusion, we recover all the states of Heisenberg and 
Schrodinger minimum uncertainty. Thus, the complete structure of minimum quantum 
uncertainty is derived by considering the osmotic contribution as the correct functional to 
be extremized. 
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5. Discussion and outlook 

The analysis carried out in the previous sections implies two main considerations. The 
first one is the emergence of the crucial role played by the osmotic uncertainty product in 
determining the fundamental features of the quantum mechanical uncertainty. 

The other main implication is that the variational approach lends itself naturally to a 
further generalization. Namely, it would be very interesting to exploit the dynamics directly 
as a constraint in the variational procedure to minimize the osmotic uncertainty. The 
physical motivation for such a construction is that it might allow to discuss the behaviour 
and the structure of quantum uncertainty for general non-harmonic systems. For instance, 
it would be interesting to verify whether the Heisenberg uncertainty product exhibits local 
minima depending on the choice of the potential. 

To implement the suggested generalization one should consider the functional 

rt 2 



C[p,v] = f 2 dtm 2 Aq 2 Au 2 
Ju 



rt2 r+co 

+ I dt I dx\ Xo(t)p + Xi(x, t) [d t p + d x (pv )] 

-,2 X2 



tl 



+ \ 2 {x,t) 



ft o d xVP d x V 

d t v + vd x v - t—^o x — J=- H 

2m z y/p m 



(40) 



The above defines the most general osmotic uncertainty functional in quantum mechan- 
ics, including the normalization, the kinematical and the dynamical constraints. Extrem- 
izing it as usual against smooth independent variations of p and v one obtains, as in the 
kinematically constrained scheme, two coupled equations. However, in this case the diffi- 
culty of identifying and/or eliminating the unknown Lagrange multipliers has prevented so 
far an explicit exact solution of the problem. Work is in progress, to obtain both numerical 



and analytical approximated solutions, and we will report on it elsewhere [12 . 
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